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INTRODUCTION 


The  characterization  of  impact  sensitivity  of 
explosive  materials  using,  for  example,  the  Picatinny 
Arsenal  impact  testing  machine  involves  some  decisions 
on  how  to  deal  with  data  provided  by  the  machine.  The 
raw  data  normally  consists  of  the  number  of  samples 
which  have  fired,  and  the  number  of  "no-fires"  for  a 
given  sample  size  N  at  each  drop  height.  In  order  to 
reduce  this  data  for  the  purpose  of  comparing  two  dif¬ 
fering  materials  or  to  decide  whether  or  not  a  particular 
material  has  changed  its  performance  (with  respect  to 
impact  measurements)  after  some  treatment,  it  is  neces¬ 
sary  to  choose  a  particular  statistical  method.  The 
choice  of  method  will,  in  turn,  have  an  effect  on  the 
size  of  sample  used  and  on  the  number  of  heights  at 
which  tests  are  made.  It  would  be  very  convenient  to 
have  one  standard  method  that  would  enable  easy  compari¬ 
son  of  the  results  of  different  workers,  or  of  measure¬ 
ments  made  at  different  times  and  places;  however, 
this  is  not  possible.  Test  procedures  which  are  appro¬ 
priate  to  the  detailed  comparison  of  only  two  materials 
may  not  be  suitable  where  it  is  necessary  to  character¬ 
ize  samples  with  respect  to  a  large  number  of  variables, 
or  where  limited  quantities  of  material  are  available. 

The  need  for  statistical  methods  in  the  reduction 
of  impact  data  comes  not  only  from  a  general  lack  of 
reproducibility  of  testing  techniques,  but  from  the 
following  considerations:  If  one  considers  a  number 
of  tests  on  a  sample  at  one  drop  height,  the  question 
occurs  as  to  why  one  drop  results  in  a  fire  and  another 
in  a  no-fire,  since  one  is  supposedly  dealing  with  the 
same  material  and  reasonable  precautions  have  been 
taken  to  make  the  test  conditions  the  same.  The  same 
question  arises  when  one  finds  that  a  test  on  one 
sample  results  in  a  fire  at  one  drop  height  and  a 
second  test  results  in  a  no-fire  at  a  greater  drop 
height  (where  the  energy  input  to  the  sample  is  un¬ 
doubtedly  higher) .  One  can  only  conclude  that  indivi¬ 
dual  samples  of  the  material  are  not  in  fact  the  same, 
but  differ  in  some  respect  which  will  result  in  differ¬ 
ing  impact  behavior. 
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This  report  is  not  concerned  with  the  physical 
or  chemical  reasons  for  these  differences  (one  plaus¬ 
ible  physical  model  is  given,  for  example,  in  the  paper 
by  Levy  et  al.  Reference  1).  It  is  merely  accepted 
that  there  exist,  within  a  sample  of  any  given  size, 
variations  in  sensitivity.  Figure  1  shows  a  distribu¬ 
tion  of  sensitivities  for  a  particular  material  such 
that  a  test  on  a  sample  at  a  drop  height  such  as  IhJ 
in  Figure  1  would  result  in  a  fire  if  the  sample  be¬ 
longed  to  the  part  of  the  distribution  to  the  left  of 
the  vertical  line  at  h  (i.e.,  if  its  critical  height 
were  less  than  the  test  drop  height)  and  result  in  a 
no-fire  if  it  belonged  to  the  larger  part  of  the  dis¬ 
tribution  to  the  right  of  that  line.  Tests  on  a  suf¬ 
ficiently  large  number  of  samples  at  a  large  number  of 
heights  would  result  in  the  typical  cumulative  curve 
shown  in  Figure  2. 

The  methods  given  below  for  the  analysis  of  impact 
data,  whether  they  involve  trials  at  only  one  or  at 
more  than  one  height,  or  whether  the  heights  are  chosen 
depending  on  the  results  of  the  testing  as  it  proceeds 
(as  in  the  Bruceton  test) ,  are  eventually  designed  to 
give  information  about  the  distribution  of  critical 
heights.  The  demonstration  of  differences  between 
materials  or  of  changes  in  a  material  as  a  result  of 
some  treatment  depends  on  showing  that  there  are  dif¬ 
ferences  between  critical  height  distributions. 

In  order  to  categorize  different  materials  it  is 
possible  to  use  properties  of  the  critical  height  dis¬ 
tribution  quite  arbitrarily:  If  trials  are  made  at 
only  one  height,  th**  estimate  of  the  percentage  of  the 
distribution  of  critical  heights  which  lies  at  heights 
smaller  that  the  test  height  might  be  used  as  some 
measure  of  sensitivity.  If  more  is  known  about  the 
distribution  from  trials  at  a  variety  of  heights,  the 
mean  of  the  distribution  might  be  used  as  a  measure  of 
"sensitivity"  and  some  measure  of  the  width  of  the 
distribution  (e.g.,  its  variance)  as  a  measure  of 
"predictability"  of  the  samples. 


In  the  discussion  of  methods  that  follows,  these 
assumptions  have  been  made: 

1„  It  is  assumed  that  the  main  purpose  of  the 
testing  is  to  detect  possible  differences  between 
materials,  and  that  the  categorizing  of  materials  with 
regard  to  their  relative  impact  sensitivities  is  of 
only  secondary  importance. 

2.  It  is  assumed  that  errors  or  uncertainties 
in  the  testing  procedure  will  have  less  effect  on  the 
uncertainties  in  the  results  than  those  arising  from 
small  sample  sizes  and  from  the  generally  broad  cri¬ 
tical  height  distributions. 

3.  It  is  sometimes  assumed  that  the  important 
parameter  in  impact  testing  is  energy  input  to  the 
sample,  and  that  drop  heights  may  be  changed  without 
affecting  the  results,  provided  that  the  total  potential 
energy  remains  the  same.  With  that  assumption,  one 
would  use  critical  energy  distribution  and  not  critical 
height  distribution.  The  assumption  is  made  here  that 

a  constant  hammer  mass  is  used  and  that  drop  height  is 
the  only  variable. 

Some  of  the  more  important  methods  will  be  given 
below.  They  are  standard  methods  and  are  discussed  in 
derail  elsewhere.  References  2  and  3  are  good  sources 
for  information  and  further  references.  Thus  only  out¬ 
lines  of  the  methods  and  some  of  their  advantages  and 
disadvantages  will  be  given. 


TESTS  AT  A  SINGLE  DROP  HEIGHT 

If  comparisons  between  materials  are  to  be  made 
from  tests  con^’^ted  at  a  single  drop  height,  the 
simplest  procedure  is  the  following: 

1.  A  height  must  be  selected  which  falls  at 
some  place  of  interest  in  the  critical  height  distribu¬ 
tion.  The  test  height  is  usually  selected  from  previous 
measurement. 

2.  A  sample  of  ft  trials  is  run  which  result  in 
Cn]  fires  and  [N  -  nj  no-fires. 
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3.  The  percentage  fire,  p  =  100  xn/N,  is 
calculated.  The  number  p  gives  an  estimate  of  the  per¬ 
centage  of  area  under  the  critical  height  distribution 
(see  Fig  1)  which  is  to  the  left  of  the  vertical  line 
at  the  drop  height  h. 

4.  A  standard  deviation  for  the  number  of 
fires  may  be  calculated: 


A  value  for  the  number  of  fires  obtained  for  some 
other  material  may  be  compared  with  this  value  of  n 
by  noting  that  percentage  areas  under  a  normal  curve 
for  numbers  of  fires  which  differ  by  s,  2s,  and  3s  from 
the  value  for  n  are  respectively  34%,  14%,  and  2%. 

A  more  elaborate  procedure  may  be  adopted  where  it 
is  assumed  that  the  procedure  for  getting  the  control 
sample  or  some  other  variable  such  as  time  of  prepara¬ 
tion  may  affect  results.  An  example  would  be  a  case 
where  materials  are  prepared,  tests  ere  run  with  a  fixed 
sample  size,  and  it  is  wished  to  make  the  control  sample 
a  composite  of  a  large  number  N  of  batches.  In  this 
case  the  number  of  fires  would  be  calculated  for  each 
of  the  batches  and  an  estimate  of  the  standard  devia¬ 
tion  of  n  calculated  from 


where  n  is  the  mean  of  the  values  of  n. 

Here  n  and  s  may  be  used  to  determine  a  confidence 
interval  for  the  value  of  the  mean.  The  interval  is 
obtained  from 

S  ±  t  s//nT 

where  t  is  found  from  tables  of  the  percentiles  of  the 
student  t  distribution  for  the  desired  level  of  confi¬ 
dence  and  N  -  1  degrees  of  freedom. 
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Statements  that  may  bo  made  when  comparing  results 
from  two  separate  tests,  A  and  B  (conducted,  of  course, 
at  the  same  drop  heights) ,  are,  for  example,  "The  mean 
percentage  fire  for  sample  A  does  (or  does  not)  lie 
within  one  (or  more)  estimated  standard  deviations  for 
sample  B,"  or:  "The  mean  of  sample  A  lies  (or  does  not 
lie)  within  some  specified  confidence  limits  for  the 
mean  of  sample  B." 

Tests  at  a  single  drop  height  may  be  made  if  it  is 
assumed  that  a  height  is  selected  at  which  changes  in 
the  critical  height  distribution  will  result  in  a 
change  in  the  percentage  fire.  In  order  to  examine 
more  of  the  critical  height  distribution  it  is  neces¬ 
sary  to  make  measurements  at  more  than  one  height. 


MEASUREMENTS  AT  MORE  THAN  ONE  HEIGHT 
Karber  Method 


This  method  assumes  that,  for  one  material,  meas¬ 
urements  of  percentage  fire  have  been  made  at  a  number 
cf  different  heights,  and  yield  results  in  the  form  of 
an  estimated  mean  critical  height  and  an  estimated 
standard  deviation.  On  the  assumption  that  the  critical 
height  distribution  is  normal  (contrary  to  experiment) , 
values  for  the  mean  critical  height  and  standard  devia¬ 
tion  completely  specify  the  performance  of  the  sample 
in  drop  testing.  The  mean  may  be  used  as  some  measure 
of  the  sensitivity  and  the  standard  deviation  as  some 
measure  of  the  predictability  of  impact  behavior  of  the 
sample.  In  cases  where  materials  are  assumed  to  be 
different,  the  Karber  method  is  useful  in  ordering 
materials  as  regards  sensitivity.  It  is  not  useful  in 
detecting  small  differences  between  materials. 

Briefly,  the  calculation  assumes  that  the  differ¬ 
ences  in  percentage  fire  at  two  heights  (for  example, 
h  and  h1  in  Figure  2)  gives  the  percentage  of  the 
critical  distribution  lying  between  the  class  limits  b 
and  h1,  and  that  sufficiently  large  numbers  of  tests 
are  conducted  at  each  height  so  that  the  variance  of 
percent  fire  at  that  height  may  be  neglected.  This 
calculation  is  given  in  detail  in  Reference  2. 


Probit  Analysis 

This  method  is  discussed  in  detail  in  Reference  2 , 
Measurements  made  at  a  number  of  different  heights  are 
tested  against  an  assumed  functional  form  for  the  crit¬ 
ical  height  distribution.  While  the  method  may  be  use¬ 
ful  for  testing  the  validity  of  particular  functional 
forms  for  the  distribution  of  critical  heights,  it  is 
not  considered  to  be  of  much  use  for  comparing  materials 
or  for  detecting  differences  between  them.  Attempts  at 
establishing  theoretical  shapes  for  impact  test  curves 
(see,  for  example.  Ref  1)  have  shown  it  is  not  possible 
to  give  them  simple  functional  forms. 

"Chi-squared**  or  Goodness  of  Fit  Test 

This  test  is  discussed  in  References  2  and  3. 

Here,  comparison  is  made  between  two  materials  which 
have  been  tested  at  a  number  Njj  of  drop  heights.  At 
each  drop  height  the  samples  are  tested  for  differences 
in  the  proportion  of  "fires'*  and  "no-fires"  anc  at  each 
height  a  value  of  x2  is  calculated  from 
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"2  («'i 
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nj  +  n 
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ni  - 
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where  nx  and  n2  are  the  numbers  of  fires  for  the  two 
samples,  and  Nj  and  N2  are  the  numbers  of  trials.  (In 
some  cases  it  is  necessary  to  group  the  results  from 
heights  which  lie  near  the  extremes  of  the  critical 
height  distribution.  This  is  discussed  in  Reference  3.) 
The  resulting  values  of  x2  are  summed  for  all  heights 
Ex2*  For  a  chosen  level  of  confidence  and  the  number 
of  degrees  of  freedom  corresponding  to  the  number  of 
heights  a  value  of  x^  _  a  is  obtained.  A  value  of 

Ex2  larger  than  xf  _  a  indicates  that  the  samples  can¬ 
not  be  said  to  be  the  same  in  regard  to  their  propor¬ 
tions  of  "fires"  and  "no-fires." 
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This  test  has  the  advantage  that  specimens  may  be 
compared  without  making  any  assumptions  about  the  dis¬ 
tribution  of  critical  heights.  If  the  drop  heights 
cover  a  sufficiently  wide  range,  the  test  shou  1  be 
sensitive  to  changes  in  that  distribution.  One  does 
not  obtain,  however,  any  information  about  the  distribu¬ 
tion,  or  any  means  of  placing  samples  in  order  of  sen¬ 
sitivity. 


Bruce ton  or 


and  Down"  Methods 


These  methods  are  also  described  in  Reference  1 
(Sec.  10.6).  They  enable  one  to  make  estimates  of  the 
mean  and  standard  deviation  of  the  critical  level  dis¬ 
tribution.  Samples  are  tested  and  the  drop  height  is 
adjusted  for  each  succeeding  test,  depending  on  the 
result  of  the  last.  These  tests  have  the  advantage 
that  they  require  no  knowledge,  or  at  least  only  a 
very  rough  estimate,  of  the  mean  of  the  critical  height 
distribution.  However,  like  the  Karber  method,  they 
assume  a  normal  distribution,  and  provide  no  simple 
way  of  testing  two  samples  for  differences. 


SUMMARY 

If  one  sets  out  to  detect  minor  differences  in 
impact  sensitivity  between  materials  or  changes  in  a 
material  the  y2  or  "goodness  of  fit"  test  recommends 
itself.  It  does  not  involve  an  assumption  about  the 
functional  fcsm  of  the  critical  height  distribution. 

In  addition,  unlike  the  Bruceton  methods  or  tests  in¬ 
volving  single  drop  heights,  it  is  sensitive  to  dif¬ 
ferences  in  any  part  of  the  distribution.  However, 
where  a  large  amount  of  testing  or  limited  amounts  of 
material  are  involved,  such  methods  may  have  to  be  used. 

The  main  disadvantage  of  the  y2  method,  that  it 
fails  to  provide  any  means  of  categorizing  materials 
as  regards  relative  sensitivity,  may  be  ■'vercome  by 
combining  it  with  one  of  the  simpler  tests  such  as  the 
Karber  method.  The  x2  and  Karber  tests  have  been  coded 
in  FORTRAN  IV  for  the  C.D.C.  6600  computer.  A  listing 
of  the  program  is  available  from  the  Explosives  Labora¬ 
tory  on  request. 
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LLi .  LM  4-1  UX 


The  general  statistical  techniques  discussed  here 
can  be  applied  to  other  sensitivity  tests.  For  example, 
the  so-called  "arbor  test"  can  be  treated  by  the  same 
statistical  methods.  In  the  "arbor  te3t"  an  uniaxial 
load  is  slowly  applied  to  the  explosive  mixture  in  a 
suitable  container  geometry  until  explosion  occurs, 
with  the  force  applied  noted  at  time  of  explosion. 
Instead  of  a  distribution  of  critical  impact  drop 
heights  for  a  given  explosive  composition  one  obtains 
a  distribution  of  critical  applied  forces  for  the  ex¬ 
plosive  composition.  Once  again  the  mean  critical 
force  and  standard  deviation  about  this  mean  indicate 
relative  sensitivity  and  predictability  of  the  explo¬ 
sive  composition.  How  closely  the  distribution  of 
critical  forces  approaches  a  normal  distribution  again 
will  help  evaluate  the  relative  merit  of  the  various 
statistical  techniques  discussed  above. 
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Cumulative  "percent  fire"  curve  (solid  line)  obtained  from 
on  large  number  of  samples  at  various  heights.  The  dashed 
sents  the  corresponding  distribution  of  critical  heights. 
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